A new improved optimization of perturbation theory: applications to 
the oscillator energy levels and Bose-Einstein critical temperature 
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Improving perturbation theory via a variational optimization has generally produced in higher 
orders an embarrassingly large set of solutions, most of them unphysical (complex). We introduce 
an extension of the optimized perturbation method which leads to a drastic reduction of the number 
of acceptable solutions. The properties of this new method are studied and it is then applied to 
the calculation of relevant quantities in different (j> 4 models, such as the anharmonic oscillator en- 
ergy levels and the critical Bose-Einstein Condensation temperature shift AT C recently investigated 
by various authors. Our present estimates of AT C , incorporating the most recently available six 
and seven loop perturbative information, are in excellent agreement with all the available lattice 
numerical simulations. This represents a very substantial improvement over previous treatments. 
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The variationally improved or optimized perturbation 
based on the linear S expansion (LDE) 0, 0, 0] is a well- 
used modification of the usual perturbation theory, based 
on a reorganization of the interacting Lagrangian such 
that it depends on an arbitrary mass parameter, to be 
fixed by some optimization prescription. In D = 1 
theories, such as the quantum mechanical anharmonic 
oscillator 4], the LDE turns out to be equivalent Q 
to the "order-dependent mapping" (ODM) resummation 
method 0. At the same time, the principle of mini- 
mal sensitivity (PMS) optimization, which takes ex- 
trema with respect to the mass parameter, is equivalent 
at large orders to a rescaling of the adjustable oscilla- 
tor mass with perturbative order, which can essentially 
suppress the factorial large order behavior of ordinary 
perturbative coefficients. This appropriate rescaling of 
the adjustable mass gives a convergent series 0, @ e.g. 
for the oscillator energy levels Q and related quantities. 
Any physical quantity whose ordinary perturbative se- 
quence is available can then be evaluated to an order 5 k 
using simply modified Feynman rules as implied by the 
following formal substitution valid for a scalar field the- 



ory: 



<(i-*) 1/2 ; 



(i) 



where lo and g are the mass and coupling respectively. 
Note that for the D = 1 quantum mechanical anharmonic 
oscillator described by a g<p 4 theory no renormalization 
is needed while for D > 1 models, the parameters 
g and lo in Eq. are to be considered implicitly bare 
parameters, and the procedure can be made in this case 
fully consistent 0,0 with the renormalization program of 
ordinary perturbation theory. In particular, appropriate 
renormalization takes into account properly any (mass or 
field) anomalous dimensions when the latter are relevant. 

Now, a definite drawback of the optimization prescrip- 
tion is that it involves minimization of a polynomial equa- 
tion of order k in the relevant mass parameter lo at per- 
turbative order 5 k , such that more and more solutions, 



most of them being complex, are to be considered when 
increasing the order. This non-uniqueness of the opti- 
mized solution requires extra choice criteria, and thus 
may seriously obscure the interpretation and the con- 
vergence towards the correct result in many non-trivial 
cases where the exact non-perturbative result is totally 
unknown. Moreover, the fact that most solutions are 
complex is embarrassing, as one has to invoke still an 
extra criterion to select a (supposedly correct) real re- 
sult. For some of the simplest models where the method 
applies, like the oscillator energy levels typically, fortu- 
nately all of the complex optimization solutions have ac- 
tually small imaginary parts and rapidly decreasing as 
the perturbative order increases (see e.g. Bcllct et al 
in refs. |||), so that the convergence properties are not 
much obscured by this inconvenience of the PMS. But 
in less trivial situations, the imaginary parts of the PMS 
solutions may be large (see e.g [Sj [l(J, ) and thus their 
physical intepretation unclear. 

In this paper, we propose a simple generalization of 
the PMS criterion as performed on the LDE series, which 
turns out to lead to a drastic reduction of physically ac- 
ceptable real optimization solutions at each successive 
perturbative order in all the physical cases we have ap- 
plied it to. First, we treat the oscillator ground state 
energy, both in the l&ige-N case (for the vector O(N)- 
symmetric </> 4 model) and the ordinary oscillator (scalar, 
N = 1) <fi 4 model. We then apply it to a less trivial 
and more interesting problem associated with the break- 
down of perturbation theory near a critical point, namely 
the evaluation of the critical transition temperature for 
a dilute, weakly interacting homogeneous Bose gas. This 
has been the source of controversy for many years and 
recently several independent groups have provided com- 
parable evaluations of the critical temperature. The rele- 
vant field theoretic framework is a </> 4 D = 3 model (after 
dimensional reduction) with an 0(2) symmetry (see e.g. 
Refs.fT^. liH ITl. IT5I Hr| for reviews). Here, for com- 
pleteness we also consider the large- TV limit of the O(N) 
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symmetric model case where the exact next-to-leading 
1/N result is known analytically. In all these cases, the 
new method seems to give excellent approximations in 
comparison with the standard PMS ones. 

I. BASIC METHOD AND THE OSCILLATOR 
ENERGY LEVELS 

Let us start with the basic perturbative series of the 
oscillator ground state energy level as described by a gcj) 4 
D = 1 model with mass lj: 

n 
9=1 

with ci = 3/4, C2 = 21/8, etc. . . 0. As discussed in the 
introduction, the LDE procedure is implemented by the 
substitutions 

Lo^Loil-5) 1 ' 2 , g^gS, (3) 

into the perturbative series Eq. JSJ) then reexpanding 
the latter to order k in the new expansion parameter 
5. Next, 5 is set to the value 6=1, such as to re- 
cover the original (massless) theory, while at any finite 
order k there remains a dependence on to in the LDE 
result E^(lj,S = 1). The standard, mostly used opti- 
mization criterion is the principle of minimal sensitivity 
(PMS), requiring at each successive perturbative orders 
k: dE^ k \tu,8 = l)/duj = for optimal ui values. The 
modification (generalization) that we propose here is first 
to introduce extra arbitrary parameters, starting at order 
two with one more parameter: 

LO^Loil-Sy^il + il-^S} 1 / 2 , (4) 

such that the modified Lagrangian still interpolates be- 
tween the free field (massive) theory for 6 — and 
the original (massless) theory for 6 = 1. Second, 
now the PMS criterion is generalized by requiring both 
dE^/du) = 0, and d 2 E^/duj 2 = 0, which gives a sys- 
tem of two equations to be solved simultaneously for a 
and lj: For a = 1, E^ 2 \lo) has no real minimum, but an 
inflection point with an almost horizontal tangent, and a 
small value of a — 1 makes this tangent horizontal, remov- 
ing the embarrassment of complex extrema. At higher 
orders, the generalization is easily done with additional 
parameters and additional vanishing of higher derivatives 
of E^ (lj). For example, at third order, we introduce two 
parameters (a and 6) as 

Lj^Lj(l-6)^ 2 [l + (l-a)S + bS 2 } 1/2 , (5) 

together with the extra requirement on the third deriva- 
tive: d 3 E^/du 3 = 1 . 



1 For completeness, note that we generalize Eq. J^J at arbitrary 
higher orders by simply adding to the expression within the right- 



The results of the new method for the (scalar case) 
oscillator energy level are shown in Table I, up to order 
nine, where higher order generalizations of Eq. 10) are 
still reasonably tractable upon using efficient poly nomial 
equation solver (we used mainly Mathematica |l7j1 2 . The 
real solution is unique (at least up to the highest pertur- 
bative orders that we analyzed, as given in Table I) and 
we also indicate for later discussion the corresponding 
values obtained for the extra interpolation parameter a as 
defined in Eqs. |@J, (JSJ. (Note that actually the (unique) 
real optimization solution gives values for the other pa- 
rameters b, ... in e.g. Eq. J5J) at higher orders, which 
are almost always small with respect to a, more on this 
later). 

We see that after approaching the exact answer to 
within 4 x 10~ 3 at fourth order, the approximation be- 
comes worse in a way strongly reminiscent of the behavior 
of an asymptotic expansion with alternating signs eval- 
uated at a finite value of its parameter. This is clearly 
related to the fact that the original series in Eq. has 
factorially growing coefficients at large orders, and thus 
our alternative method apparently looses the property 
of compensating this divergent behavior by an appro- 
priate rescalingja, |(| of the optimized mass parameter, 
unlike the standard LDE-PMS. More precisely, one can 
see from the third column of Table I that the values of 
the optimized mass parameter lj, after increasing regu- 
larly with the perturbative order for the first few orders, 
appear rather to have a slower increase for the highest 
perturbative orders we could consider. 

Nevertheless, we can try to better exploit these results 
by standard resummation methods applied on the ob- 
tained perturbative sequence, typically considering sim- 
ple Pade approximants. Accordingly, let us define the 
following perturbative series: 

£b) =El + - E <i) x<1 > ( 6 ) 

9=1 

where E q denotes the results for the ground-state en- 
ergy at orders q in Table I, so that Eq. © simply gives 
E(p) = E p f or x = 1. Now a simple (3,3) Pade approx- 
imant of Eq. © using the first seven orders of Table I 
gives 0.420841, and a (4,4) Pade using all the numbers 
of Table I gives 0.420838, still an improvement, within 
8 x 10 -5 of the exact answer, in spite of the superficially 
disastrous results of the eighth and nineth orders. 

Next, the same analysis is performed for the large- iV 
approximation of the vector 0(N) symmetric oscillator, 
which consists of considering only the "cactus" Feyn- 
man graphs. As is well-known, the resulting perturba- 



handed bracket of Eq. a term bj c 5 k at each order k, where 
6fe are new parameters. 
2 It is clearly becoming more cumbersome and CPU consuming to 
apply such a method at very large perturbative orders due to the 
increase in the number of parameters at successive orders. 
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TABLE I: Improved-PMS results for the oscillator ground 
state energy level, at different orders k, with the correspond- 
ing values of the optimal mass parameter u> and the main 
interpolation parameter a. E exact = 0.42080497 . . . (4g) 1/3 . 



k 


Eipms/(4 9 ) 1/3 


Ul 


a 


1 


0.429268 


1.82 




2 


0.418483 


2.04 


1.05 


3 


0.422341 


2.18 


1.08 


4 


0.419138 


2.28 


1.09 


5 


0.423496 


2.37 


1.11 


6 


0.41523 


2.45 


1.12 


7 


0.436015 


2.51 


1.15 


8 


0.380812 


2.59 


1.15 


9 


0.647259 


2.60 


1.25 



TABLE II: Improved-PMS results for the oscillator ground 
state energy level at large N, for different orders k, 
with the corresponding values of the optimal mass pa- 
rameter ui and interpolation parameter a. E exac t = 
0.4292678409... (4 g/3) 1/s . 



k 


m l ( a /o\l/3 

EiPMS/ykg/S) ' 


Ul 


a 


1 


0.429268 


1.26 




2 


0.429589 


1.37 


0.96 


3 


0.429400 


1.43 


0.92 


4 


0.429326 


1.47 


0.89 


5 


0.429296 


1.50 


0.87 


6 


0.429282 


1.52 


0.85 


7 


0.429275 


1.54 


0.84 


8 


0.429272 


1.55 


0.83 


9 


0.429270 


1.56 


0.82 



tive series for the ground state energy of a form similar 
to Eq. © can be obtained simply in this case to arbi- 
trary orders by solving the large- TV gap equation exactly, 
and expanding in a perturbative series in g the corre- 
sponding exact large-A^ expression of the ground state 
energy. When written with the same normalization as in 
Eq. |2l the first few coefficients of this series are c\ = 1/4, 
c-2 = 1/4, C3 = 1/2,.... However, an important differ- 
ence with the scalar case oscillator is that the resulting 
perturbative series has no factorially growing coefficients 
at large orders and a finite convergence radius. As one 
can see, the method performs very well in this case, ap- 
proaching rapidly very close to the exact solution when 
the order increases. Again, there is only one real pos- 
itive solution at each order. One can also notice from 
the last column of Table II the more regular behavior 
of the extra interpolation parameter a, as compared to 
Table I, which now appears to tend towards a constant 
value. (The other parameters b, etc. . . corresponding to 
the analysis in Table II are numerically quite negligible 
with respect to a, so that we do not give their explicit 
values) . 

In order to better understand why the method seems 
not so efficient when applied to a factorially divergent se- 
ries, like Eq. @ for the standard (scalar case) oscillator, 
let us consider now the even simpler functional integral 
of euclidian (f> theory in zero dimension: 

I{g, m) = J dx exp{-m 2 x 2 /2 - g x 4 /A} , (7) 

As is well-known, this can be expanded in a perturba- 
tive series in g with alternating signs, factorially growing 
coefficients at large orders: 

(8) 



Now the exact result of the integral for m = is 

1(9, 0) = 2-1 ^^9- 1/4 - 1.0227656721 . . .g-V* , 

which we can compare with the results of the applica- 
tion of the method: 0.954325, 1.0865771, 0.94213824, 
1.244391803, 0.7038846 respectively for orders one to five. 
Then, there is no real solution at order six. Nevertheless, 
a Pade approximant f[i,i] (thus constructed from the 
first three orders) still gives 1.01754. Thus, the behavior 
is quite similar to the one of the (scalar) oscillator where 
there are also no physically acceptable solutions at higher 
orders. (However, the quality of the results from Pade ap- 
proximants of higher orders appears to deteriorate more 
rapidly than in the previous oscillator case: For exam- 
ple, a P\2,2] thus constructed from the first five orders 
gives 1.05677). This anyway confirms that the general- 
ized method is not appropriate to turn a factorially di- 
vergent large order behavior into a more convergent one, 
while it definitely improves the LDE convergence when 
starting from a less divergent original series. 

II. CRITICAL THEORY AND BEC CRITICAL 
TEMPERATURE 

We now turn to the more recent and challenging prob- 
lem of the BEC critical temperature evaluation. We first 
recall that, to 0{a 2 n 2 / 3 ), its functional form was found 
to be [3 

T c = T {l + Cl an 1/3 + [c^ln{an 1/3 ) + 4}a 2 n 2/3 } , (10) 

where To is the ideal gas condensation temperature, a is 
the s-wave scattering length, n is the density and ci, c' 2 , d{ 
are numerical coefficients. Some early analytical predic- 
tions included the self-consistent resummation schemes 
[l4j (ci ~ 2.90), the l/N expansion at leading order 
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(ci ~ 2.33) and at next to leading order 0] (ci ~ 1.71) 
and also the LDE at second order [3 (ci ~ 3.06). The 
numerical methods include essentially lattice simulations 
(LS). The most recent LS results are reported by the au- 
thors of Ref. EH (ci = 1.29 ± 0.05) and of Ref. HJ 
(ci = 1.32 ± 0.02). Very recent analytical studies predict 
ci = 1.27 ± 0.11 HU. The problem is that these coeffi- 
cients (except c' 2 ) are sensitive to the infrared physics at 
the critical point and so, no perturbative approach can 
be used to compute them. At the critical point, one can 
describe a weakly interacting dilute homogeneous Bose 
gas by an effective action analogous to a 0(2) scalar field 
model in three-dimensions given by 

S*=J d 3 x Q(V0) 2 + \^ + ^ , (11) 

where <fi is a two-component real scalar field. The param- 
eters ui and g are related to the original parameters of the 
nonrelativistic action by to = —Imp, and g = A&iramT 
|l5l ll6T ] where \i represents the chemical potential, m 
the atomic mass and T the temperature. The leading 
order coefficient of the critical temperature shift can be 
expressed as ^3 

d = -(256/Af)7r 3 [C(3/2)]- 4 / 3 K , (12) 

where gn — A(0 2 ) = {4> 2 ) g — (<t> 2 )o- The subscripts g 
and mean that the field fluctuations are to be evalu- 
ated in the presence and in the absence of interactions 
respectively. 

The implementation of LDE within this model is re- 
viewed in previous applications 0, 0, 0, 0|. Let us 
first consider the LDE in the large- N limit extending the 
relevant 0(2) model to a 0(A)-symmetric one. The in- 
terest of such an approximation is that the exact c\ result 
is known by direct evaluation |l5j . The original pertur- 
bation series in this case can be written as 

<w~£ + £e«(-£)'. <-> 

i— 1 v / 

where the perturbative coefficients are given by 
3 / 1 Y f°° z 2 

(14) 

with 



A(z) = — arctan— , (15) 

and z — k/uj, and can be calculated anal ytic ally to ar- 
bitrary precision using, e.g., Mathematica | l7| numerical 
integration. By applying the generalized LDE procedure 
as described above in Eqs. J3J, (JSJ, etc. . . up to order 
S 9 , one obtains the results shown in Table III (ci./pms), 



TABLE III: Improved-PMS (IPMS) versus standard results 
for the large-TV BEC AT C at different orders k. ci )exa ct = 
2.32847 . . . 



k 


ClJPMS 


a 


CifiestPMS 


2 


2.163 




2.163 


3 


1.85 


1.44 


1.88 ±0.17/ 


4 


1.93 


1.69 


1.96 


5 


2.00 


1.78 


1.91 ±0.03/ 


6 


2.04 


1.83 


1.94 


7 


2.08 


1.89 


1.93 ±0.0015/ 


8 


2.10 


1.88 


1.935 


9 


2.12 


1.89 


1.95 


10 


2.14 


1.91 


1.95 ±0.03/ 



together with those obtained by the standard PMS op- 
timization procedure, the latter for the best converging 
family of solutions Ci^estPMS (see Ref. [l(j for details on 
the latter). 

As illustrated and discussed in details in [13, the con- 
vergence of the standard LDE-PMS method when ap- 
plied to the exact large- ./V series Eq. (fl3|) is actually very 
slow (leaving apart the accidentally good result at first 
order), so that only at very large perturbative orders ~ 50 
does the optimized values approach reasonably close to 
the exact large- N value, c\ ~ 2.32847... Thus the advan- 
tage of the improved PMS approach is not obvious in this 
case, since it is evidently less algebraically tractable than 
the standard method at very large perturbative orders. 
Nevertheless, one can see from the improved-PMS results 
in Table III that the new method performs much better, 
giving always a unique real solution and seemingly con- 
verging much faster towards the expected result. In fact, 
as motivated in , m order to study the eventual LDE 
convergence properties of the large- N expression of c±, 
it is more instructive to consider an approximated form 
of the corresponding large- N perturbative series which 
takes into account only the relevant infrared limit of the 
auxiliary field propagator (see Ref. ^3 f° r more details). 
This has the advantage of giving a simple geometric series 
with exact coefficients: 

(A?\W Nid gN k ^ f gNY 

i— 1 v ' 

where Gi = [(647r 2 )(87r) 1 ] -1 , such that the straightfor- 
ward resummation of (|16J) is well defined for the relevant 
limit id — > 0, which accordingly can be reached smoothly 
in contrast with the genuine large- A series Eq. (I13fl . The 
alternative method results applied on the original series 
(|l(j[l are shown in Table IV, together with those obtained 
by the standard PMS optimization procedure, the latter 
for the best converging family of solutions. 

In that case, the convergence of the standard PMS 
method is faster in comparison of Table III, but the im- 
proved PMS method performs even better: starting from 
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TABLE IV: Improved-PMS (IPMS) versus standard results 
for the large-iV BEC AT C at different orders k. c\^ xa ct = 
2.32847 . . . 



k 


ClJPMS 


a 


Cl,bestPMS 


2 


2.852 




2.852 


3 


2.36771 


1.49 


2.444 ±0.276/ 


4 


2.34451 


1.53 


2.244 ± 0.20/ 


5 


2.33686 


1.56 


2.397 ±0.079/ 


6 


2.32847 


2.0 


2.333 ± 0.08/ 


7 


2.32847 


2.0 


2.298 ± 0.06/ 


8 


2.32847 


2.0 


2.342 ± 0.04/ 


9 


2.32847 


2.0 


2.324 ±0.036/ 



order six and beyond, the exact result is always obtained 
as the unique physically acceptable solution. It is inter- 
esting to trace what happens in more details. In fact, 
the solutions found at LDE orders k > 6 by applying the 
procedure isa = 2, b = c = d = . . . = 0, which is easily 
seen by comparing with Eqs. J5| to correspond to a 
basic interpolation of the form 

lu -» w(l - 2 8 + S 2 + ) 1/2 = uj(1 - S) , (47) 

Then, applying the substitution (|17|) on any geometric se- 
ries instead of the standard LDE substitution JIJ canon- 
ical for a scalar theory, one can easily see from simple 
algebraic properties that the improved PMS solution al- 
ways reproduces the exact solution, and this at any arbi- 
trary LDE order k > 2. This case exhibits a spectacular 
improvement over the standard LDE results, which con- 
verges only rather slowly and with non-zero (albeit small) 
imaginary parts, a source of much frustration. A pos- 
teriori, there is nothing particularly remarkable in this 
result which is essentially an algebraic accident of the 
non-canonical substitution Eq. I|f 7|) followed by the LDE 
when performed on a simple geometric series. However, 
what is perhaps more interesting is that our improved- 
PMS procedure is in that way guessing a more appro- 
priate value of the rescaling power within the LDE sub- 
stitution Ansatz Eq. Alternatively we could have 
parametrized the basic LDE substitution according to 

oj^oj(l-5y , (18) 

with an arbitrary power 7 to begin with, and then look 
for the best value of 7 such that the LDE series con- 
verges faster towards the exact result, which is clearly 
the case for 7 = 1 in this large- N case. Now, it is worth 
noting that considering an arbitrary power coefficient 7 
according to Eq. I|18l) when applied to the BEC series 
turns out in practice to be essentially equivalent to mod- 
ifying the simplest LDE substitution formula Eq. by 
introducing the relevant critical exponent: 

lu -> lu(1 - 5) 1/u1 ' , (19) 



where to' = 2f2/(2 — rj), 77 is the anomalous dimension 
of the critical propagator ~ l/p 2- * 7 and fl = 0'(g c ), g c 
being the critical coupling (see e.g. Ref. O)- This 
renormalization group inspired modification Eq. I|19|) 
of the standard LDE Eq. O ) is indeed the approach 
followed e.g. in Refs. |2ll |22|. where numerical values 
of lu' as obtained by different methods (including the 
variational perturbation theory j23J) are used in Ansatz 
(|19|l prior to an (otherwise standard) PMS optimization. 
Now in our case, a major difference is that the relevant 
exponent lu' ~ 2/ a is simply guessed by the generalized 
optimization procedure, at the same time as obtaining 
as optimized solution the relevant physical quantity 
c\. More precisely, the correct exact large- N value 
lu' = 1, fi = 1 (see e.g. Ref. 0]) is clearly guessed 
by our improved PMS procedure, with lu' — 2/ a, at 
least for the infrared approximated large-iV case as 
illustrated in Table IV. Note also that the values of 
the parameter a at successive perturbative orders are 
also clearly approaching quite closely, though more 
slowly, the correct critical value when considering 
the genuine large-A^ series, as illustrated in Table III 
(noting however that in this case the other parameters 
6, etc. . . are small with respect to a but not strictly zero). 

We will consider now mainly the physically relevant 
N = 2 case, as well as the case N — 1 and N = 4 for 
comparison with other available lattice simulation and 
analytical results. We will see that many of the previous 
results for the large- N series generalize, at least qualita- 
tively, to this less trivial case. For N = 2, the quantity 
2 )g k ^ has been first evaluated, up to order (5 4 , in Ref. 
. Recently, hig her order terms have been evaluated 
by Kastening [2l| so that, to order S 6 , the perturbative 
series can be written as 

w (-=)'• < 20 > 

i=l 

where the coefficients are given by 3 K x = 3.22174xl0~ 5 , 
K 2 = 1.51792 x 1(T 6 , K 3 = 9.66512 x 1CT 8 , K 4 ~ 
7.51366 x 1CT 9 , and K 5 ~ 6.7493 x 1(T 10 . The results of 
our alternative procedure are shown in Table V, together 
with the results from applying simple Pade approximants 
similar to the ones discussed above for the oscillator, i.e. 
defining a new perturbative series similar to Eq. © but 
where the Ei are now replaced by the values of c\ in Table 
V at successive orders k. One can notice that at orders 



3 Throughout this paper wc use the perturbative loop coefficients 
results of Ref. L'll to evaluate all these coefficients, since in par- 
ticular the latter are obtained as much as possible from exact 
integrals. There are therefore some very small differences in the 
lowest orde rs numerical results with respect to some previous 
analysis 111 1 1(1 . since the latter used numerical integration, not 
always very precise. Not e als o a trivial difference of normaliza- 
tion with respect to Ref. |21| in our defining series, Eq. 1201 . 
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k = 2, 3 and 4 the IPMS results quickly oscillate around 
the MC results. From order S 4 onwards the oscillation 
is reduced drastically and the IPMS procedure generates 
stable results which agree remarkably well with the lat- 
tice results. Furthermore, we also indicated in Table V 
the corresponding values obtained at each order for the 
parameter a, whereas as already indicated, in the gener- 
alized LDE substitutions Eq. J5J), etc. . . at higher orders 
the remaining parameters b, etc are numerically smaller. 
Accordingly, to first order approximation, 

lu -> u[l - aS - (I - a)8 2 + bS 2 (l - S) + •■ -] 1 ' 2 
~ Lu{l-8) a ' 2 , (21) 

so that in a rough approximation one has lu' ~ 2/ a. Tak- 
ing the latter approximate relation at face value would 
give e.g. for the successive orders considered in Table V: 
lu' ~ 1.03, 0.71, 0.73, 0.71, which compare not badly with 
the reported values of lu' = 0.8 ± 0.04 [H obtained 
numerically by other methods. Note, however, that the 
present method does not pretend at this stage to accu- 
rately predict in that way the relevant critical exponent 
lu', etc. Indeed, in this N — 2 case the other parame- 
ters b, etc. . . are small with respect to a but not strictly 
negligible. Nevertheless, though the essential motivation 
of our approach introducing more interpolation param- 
eters in Eq. J5J is to get rid of the unwanted complex 
optimization solutions, it is interestingly connected with 
the parametrization of corrections to scaling, the lead- 
ing correction bein g co rrectly described according to Eq. 
(|T§)) (see e.g. Ref. |2ll |2^| for further motivation of Eq. 
mU). However, to be complete one should note that, in 
contrast with the previous oscillator and large- N BEC 
series, the real IPMS solutions in Table V are not al- 
ways unique: actually, at orders k = 4 and k — 6 there 
appear a second real (but negative) solution. But these 
extra solutions can be immediately eliminated as they 
correspond to largely unreasonable values of the inter- 
polation parameters: for instance for k = 4 the extra 
real solution is: a ~ —13.9, b ~ 55.5, and Q ~ —0.02 
which gives c\ ~ —55.4, while we can expect consistent 
values of these parameters to be reasonably close to their 
corresponding large-V values, i.e. typically a should be 
not too far from its large- N value a = 2, as also sup- 
ported from the results of the above mentioned leading 
corrections to scaling analysis. 

At the same time, the results from the standard LDE- 
PMS optimization method shown in Table V, also oscil- 
late somehow around the MC result, as anticipated in 
Ref. [Tlj, but the convergence is not at all obvious: after 
approaching rather closely the lattice result at order S , 
the results of the next two higher orders depart sensibly 
from the lattice one. In the absence of any indication 
on the higher orders, we can hardly speculate but, by 
comparison with the results of Ref. 10] for the large- 
TV case, as summarized in Tables III and IV, it appears 
possible that the standard LDE ultimately converges to- 
wards the same result, but in a much slower way than 
our improved method given in the left column of Table 



V. Moreover, the latter standard PMS results will look 
much less attractive if one recalls that they have been se- 
lected among several possible complex results with large 
imaginary parts |10|. requiring an extra selection crite- 
ria 4 , in contrast to the improved version. 

One can also note the remarkable results of the Pade 
approximants based on the improved PMS results at dif- 
ferent orders, even for the lowest order one f[i,i], though 
it only uses the second and third order IPMS results. 
This is certainly not coincidental and should be mainly 
attributed to the oscillatory property of the IPMS re- 
sults. Indeed, the very same Pade approximants using 
instead the (real parts of) the standard PMS results of 
the rightmost column of Table V are far away from any 
reasonable result: this would give e.g. fri,i] = 4.26 and 

P[2,2] = 3 - 66 ' 



TABLE V: Improved-PMS {IPMS) versus standard (PMS) 
results for c\ at N = 2 in the BEC case at different orders 
k. Also shown are the Pade approximants (PA) resummation 
results applied to the improved-PMS cases. The lattice results 
are ci = 1.29 ± 0.05, 1.32 ± 0.02. 



k 


ClJPMS 


a 


PA(d t ip M s) 


&L,PMS 


2 


3.06 






3.06 


3 


0.98 


1.95 




2.45 ± 1.66J 


4 


1.426 


2.81 


P[i,i] = L347 


1.53 ±2.32/ 


5 


1.247 


2.75 


P [li2] = 1.283 ; P [2 ,i] = 1.298 


0.76 ± 2.53/ 


6 


1.300 


2.83 


P [2 ,2] = L286 


2.40 ± 1.69/ 



Next, we consider the completely similar calculation 
of the coefficient c\ as defined in Eq. I|12(l but for the 
cases N = 1 and N = 4, respectively, for which lat- 
tice calculations of ci have been recently performed in 
Ref. 26] . The original perturbative series reads like 
Eq. (|20|l with now the relevant coefficients Ki given 
by K x = 1.208097 x 10" 5 , K 2 = 5.1229853 x 10~ 7 , 
K 3 = 2.9645225 x 10" 8 , K 4 ~ 2.1084956 x 10~ 9 , and 
K 5 ~ 1.7420267 x lO" 10 , and K x = 9.6647794 x 10~ 5 , 
K 2 = 5.4645177 x 10~ 6 , K 3 = 4.09584576 x 10~ 7 , 
K 4 ~ 3.70261198 x 10~ 8 , and K 5 ~ 3.8333358 x 10" 9 , 
respectively for N = 1 and N = 4. 

Our results are given in Tables VI and VII respectively 
for N = 1 and N = 4. As one can see, the results from 
our alternative method are again showing very good con- 
vergence properties and approaching very closely the lat- 
tice results. The Pade approximants are also in excellent 
agreement with the latter. In addition, the correspond- 
ing values obtained for the parameter a, as related in first 



4 More precisely, we selected I lfl among the different LDE-PMS 
real or complex solutions the ones having the smallest /?e[il;] > 
value of the optimized mass, motivated by the fact that the exact 
solution would correspond to w — » 0. 
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approximation to the critical exponent u>' in Eq. (|18fl , ap- 
pear qualitatively very consistent to known results for to' 
for N = 1 and N = 4 HI El- On the other hand, the 
results from the standard LDE-PMS method show a be- 
havior very similar to the case N = 2: after approaching 
not too far from the lattice result at fourth order (but 
with very large imaginary parts) , the results from higher 
orders are not satisfactory with trends very similar to 
the case N = 2. 



TABLE VI: Same as Table V (improved versus standard) 
PMS results but for N = 1 at different orders k. The lat- 
tice result is 1.09 ± 0.09. 



A: 


ClJPMS 


a 


PA(ci,ipms) 


Cl,PMS 


2 


2.65 






2.65 


3 


0.817 


1.95 




2.12 ± 1.47J 


4 


1.237 


2.88 


P[i,i] = 1-159 


1.31 ± 2.06/ 


5 


1.047 


2.77 


P [li2] = 1.086; P [2A] = 1.106 


0.62 ± 2.29/ 


6 


1.114 


2.90 


P M = 1.095 


2.11 ± 1.54/ 



TABLE VII: Same as Table V (improved versus standard) 
PMS results but for N = 4 at different orders k. The lattice 
result is 1.59 ±0.10. 



k 


ClJPMS 


a 




Cl,PMS 


2 


3.75 






3.75 


3 


1.222 


1.94 




2.99 ± 1.99/ 


4 


1.665 


2.68 


P[i,i] = L589 


1.90 ± 2.74/ 


5 


1.524 


2.66 


P [li21 = 1.550; P [2ill = 1.558 


0.99 ±3.01/ 


6 


1.556 


2.74 


P [2 ,2] = 1-549 


2.90 ± 1.98/ 



We finally apply our alternative method to the other 
non-perturbative relevant coefficient c 2 in the defining 
Eq. I|1UI) . This coefficient, which is not generally con- 
sidered by most authors working on the BEC AT C prob- 
lem, was first evaluated with lattice simulations. To our 
knowledge, its sole analytical evaluation made use of the 
standard LDE-PMS U OH up to order 5 4 (five loops). 
This <3(a 2 n 2 / 3 ) coefficient appearing in the AT C expan- 
sion, Eq. I|10|l. can be written as 



4 



-5/3 6 // + I[£(3/2)]- 8 / 3 (1927r 3 K ) 2 



where b 2 



-g[C(3/2)] ■ U, 

^C(l/2)[C(3/2)]- 5 / 3 lnC(3/2) 



(22) 



is 



32tt{ [i ln(1287r 3 ) ± \ - 72tt 2 K - 96tt 2 k] 



xC(l/2)±^-X 2 -^|[C(l/2)] 2 }, (23) 

with K 2 = -0.13508335373 and g 2 K = r c where r c = 
— S(0) (Hugenholtz-Pines theorem). As before, the quan- 
tity k is easily obtained from Eq. (|20|) whereas TZ can 



be obtained directly from the perturbative evaluation of 
E(0). We refer the interested reader to Ref. |!j for the 
details and subtleties associated with this type of eval- 
uation. Thanks to the recent availability of improved 
six- loop results [24J one obtains , at N = 2, 



= -S^(0) = f +g 2 A 2 

D7T 



± w 



E 

i=3 



-<J 



Ai 



0.59775 

(24) 



where M is an arbitrary mass scale introduced by di- 
mensional regularization and A 2 = 1.40724 x 10~ 3 , 
A 3 = 8.50888 x 10" 5 , A A = 3.57259 x 10" 6 and A 5 = 
2.25332 x 10~ 7 . As discussed in Ref. it is interest- 
ing to note that the optimized u> defined by the standard 
PMS, or as well by the IPMS introduced in the present 
work, are scale independent to any order in 8 so that u 
is only g dependent. Note also that, contrary to A(0 2 ), 
r c is a divergent quantity. In this case, the optimiza- 
tion procedure must be implemented after renormaliza- 
tion as advocated in Ref. p|. Finally, let us point out 
that the potential technical difficulty associated with the 
evaluation of r c stems from the fact that, while A(cf) 2 ) 
depends on the (finite) difference S(p) — E(0), the for- 
mer de pen ds on the (divergent) S(0) only. Here, we had 
access |2j| to values of S(0) up to order S 5 (six loops) 
only. Unfortunately, at order 5 e only the results for the 
joint S(p) — S(0) contribution is known [24| . 
The results of the alternative method are shown in Ta- 
ble VIII below for M — g/3 which was the value chosen 
in the lattice evaluations H3l- As one can see, the re- 
sults of our alternative method show again a remarkable 
agreement with the numerical lattice result up to the fifth 
order (six loop) . Concerning the results of the standard 
PMS method, as shown in the right column of Table VIII, 
their behavior is again very similar to the ones for the ci 
coefficient above in Tables V-VII: after approaching very 
closely the lattice result at fourth (five-loop) order, the 
six-loop result starts to decrease a bit further. 



TABLE VIII: Improved-PMS results for c 2 ' in the N = 2 BEC 
case at different orders k. The lattice result isc 2 ' = 75.7±0.4. 



k 


C 2,IPMS 


Pade(c 2 ' JPMS ) 


C 2,bestPMS 


2 


101.2 




101.2 


3 


69.83 




94.2 ± 31.2/ 


4 


77.90 


P [M1 = 76.25 


75.0 ±41.1/ 


5 


73.46 


P M = 79.31; P i2M = 75.04 


60.4 ±41.2/ 
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III. CONCLUSION 

We have introduced a conceptually rather simple gen- 
eralization of the optimized (or variationally improved) 
perturbation method, which considerably reduces the 
number of irrelevant optimization solutions at each per- 
turbative order. We have applied this variant of the PMS 
method to the calculation of certain non-perturbative 
quantities in the simple quantum mechanical oscillator 
as well as the more challenging BEC critical tempera- 
ture determination. 

Just like the usual LDE/PMS procedures, it is a recipe 
without a formal general justification. For quantities de- 
scribed by a perturbative series which is originally facto- 
rially divergent, such as the series relevant for the oscil- 
lator, we find that this new procedure, though it helps 
in selecting only a few (if not always uniquely) among 
the many possible complex solutions, does not exhibit 
obvious convergence improvement behavior, in contrast 
with the standard method whose rigorous convergence 
is established in simple cases. Our experience from sev- 
eral toy series as well as the physical systems considered 
in this paper seems to point to the following behaviors 
which remain somewhat mysterious to us: Compared to 
the usual LDE/PMS procedures, it seems much supe- 
rior conceptually when applied to series with finite ra- 
dius of convergence or at least less divergent series, like 
the ones relevant to the BEC critical temperature. When 
one deals with an asymptotic expansion (with zero radius 



of convergence but alternating signs) typical of those ap- 
pearing in perturbative calculations for physical systems, 
it seems that the recipe breaks down at high order, but 
that at low order it gives excellent numerical results, pro- 
vided that the asymptotic expansion is not too different 
in its first few orders from a series with finite radius of 
convergence (Examples are provided by finite- N versus 
large- N theories). 

Moreover, quite interestingly the method appears in non- 
trivial cases to estimate numerically (by optimization) es- 
sentially correct values of the leading corrections to scal- 
ing behavior, since the introduction of more interpolation 
parameters in the LDE Anstaz is, to a first approxima- 
tion, equivalent to introducing the relevant critical expo- 
nent. Finally, our numerical results obtained for the BEC 
critical temperature when using the latest available per- 
turbative information are in excellent agreement with the 
numerical lattice simulations results for the all available 
finite N cases (and also for the large- N case as compared 
to the exact analytic result). 
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